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Abstract 

We  extend  our  recent  simulation  studies  where  a  screw  dislocation  in  face-centered  cubic  (fee)  Ni  was  found  to  spontaneously  attain  a 
low  energy  partially  cross-slipped  configuration  upon  intersecting  a  forest  dislocation.  Using  atomistic  (molecular  statics)  simulations 
with  embedded  atom  potentials,  we  evaluated  the  activation  barrier  for  a  dislocation  to  transform  from  fully  residing  on  the  glide  plane 
to  fully  residing  on  a  cross-slip  plane  intersecting  a  forest  dislocation  in  both  Ni  and  Cu.  The  activation  energies  were  obtained  by  deter¬ 
mining  equilibrium  configurations  (energies)  when  variable  pure  tensile  or  compressive  stresses  were  applied  along  the  [111]  direction  on 
the  partially  cross-slipped  state.  We  show  that  the  activation  energy  is  a  factor  of  2-5  lower  than  that  for  cross-slip  in  isolation  via  the 
Escaig  process.  The  cross-slip  activation  energies  obtained  at  the  intersection  in  Cu  were  in  reasonable  accord  with  the  experimentally 
determined  cross-slip  activation  energy  for  Cu.  Further,  the  activation  barrier  for  cross-slip  at  these  intersections  was  shown  to  be  lin¬ 
early  proportional  to  (d/b)[ln(\/3d/b)]1/2,  as  in  the  Escaig  process,  where  d  is  the  Shockley  partial  dislocation  spacing  and  b  is  the  Burgers 
vector  of  the  screw  dislocation.  These  results  suggest  that  cross-slip  should  be  preferentially  observed  at  selected  screw  dislocation  inter¬ 
sections  in  fee  materials. 

©  2010  Acta  Materialia  Inc.  Published  by  Elsevier  Ltd.  All  rights  reserved. 
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1.  Introduction 

There  is  increasing  recognition  of  the  need  to  incorpo¬ 
rate  physics-based  models  of  deformation  in  the  design  of 
structural  components.  While  models  for  predicting  yield 
strength  and  creep  behavior  are  beginning  to  incorporate 
significant  physics,  models  for  fatigue  and  ultimate 
strength  remain  mostly  empirical.  This  is  because  strain 
hardening  and  fatigue  resistance  are  highly  influenced  by 
dislocation  micromechanisms,  for  example  cross-slip,  and 
including  physics-based  cross-slip  processes  in  mesoscale 
simulations  (e.g.  dislocation  dynamics)  has  been  difficult. 
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The  early  work  of  Escaig  remains  the  most  widely  cited 
and  used  model  for  cross-slip  [1-4],  however,  this  model 
poses  several  difficulties  with  respect  to  quantitative  simu¬ 
lations.  The  model  is  highly  sensitive  to  choice  of  parame¬ 
ters  that  have  thus  far  been  difficult  to  obtain.  For 
example,  the  constriction  width  required  for  cross-slip  is 
unknown  but  significantly  influences  the  enthalpy  of  the 
cross-slip  process  [5],  This  difficulty  has  led  to  an  ad  hoc 
postulate  that  obstacles  always  exist  in  materials  which 
enable  sufficient  dislocation  core  constriction  under  stress, 
thereby  ensuring  cross-slip  [1],  Aside  from  being  a  pre¬ 
sumption,  this  forces  cross-slip  models  (particularly  for  sin¬ 
gle  crystals)  to  make  arbitrary  assumptions  about  the 
dislocation  obstacle  spacing.  Advances  in  atomistic  simula¬ 
tions  make  it  possible  to  gain  insights  into  the  cross-slip 
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process  and  may  serve  to  inform  mesoscale  simulations  to 
accurately  capture  the  atomic  level  physics  of  that  disloca¬ 
tion  process.  Such  an  insight  was  realized  in  our  recent 
work  [6],  where  the  possibility  of  cross-slip  at  disloca¬ 
tion-dislocation  intersections  was  selectively  examined 
using  large-scale  atomistic  simulations  that  contained  more 
than  one  dislocation.  A  glide  dislocation  that  intersects  a 
pair  of  forest  dislocations  was  modeled  and  several  condi¬ 
tions  were  shown  to  exist  where  cross-slip  can  nucleate  eas¬ 
ily.  This  provides  a  better  physical  basis  for  cross-slip  and 
the  means  to  incorporate  realistic  cross-slip  evolution  in 
larger  scale,  discrete  dislocation  dynamics  simulations. 
Further,  this  may  lead  to  a  realistic  statistical  representa¬ 
tion  of  cross-slip  during  monotonic  or  cyclic  deformation. 
Flowever,  this  prior  work  did  not  explore  the  size  and  nat¬ 
ure  of  the  activation  barrier  for  dislocation  transformation 
at  intersections. 

In  this  work,  we  show  that  a  finite  energy  barrier  exists 
for  cross-slip  at  an  intersection,  but  that  it  is  much  smaller 
than  that  obtained  from  analyses  of  the  Escaig  process. 
Molecular  statics  simulations  were  used  to  determine  the 
saddle  point  configuration  and  activation  energy  by  using 
an  applied  stress  to  change  the  extent  of  partial  dislocation 
cross-slip.  The  applied  stresses  were  compressive  or  tensile 
along  the  [111]  direction  on  the  partially  cross-slipped  dis¬ 
location  structures  obtained  in  a  prior  work  [6].  The  screw 
character  dislocation  under  study  had  a  1/2[1 1 0]  Burgers 
vector  and  the  intersecting  dislocation  had  a  1/2[I0 1]  Bur¬ 
gers  vector  on  the  (ill)  plane  and  a  120°  line  orientation 
of  [0  —1  —1]  (see  Fig.  1).  Section  2  describes  the  simulation 
technique  and  the  interatomic  potentials.  Section  3  presents 
a  review  of  the  core  structures  obtained  for  this  specific 
intersection  [6]  in  Ni,  as  well  as  core  structures  obtained 


from  atomistic  simulations  for  this  specific  intersection  in 
Cu,  the  effect  of  applied  stresses  along  the  [111]  direction 
on  the  partially  cross-slipped  core  structure,  identification 
of  the  saddle  point  configuration  and  a  determination  of 
the  energies  of  the  activated  state  for  cross-slip.  Finally, 
Sections  4  and  5  give  a  discussion  and  summary  of  the 
results,  respectively. 

2.  Simulation  technique 

The  atomistic  simulations  described  here  employed  the 
3-dimensional  (3D)  parallel  molecular  dynamics  code 
LAAMPS  [7],  developed  at  Sandia  National  Laboratory. 
A  schematic  of  the  simulation  cell  used  in  the  atomistic 
simulations  is  given  in  Fig.  1.  The  simulation  cell  was  a 
rectangular  parallelepiped  having  the  x-axis  oriented  along 
[1 1 0],  the  j-axis  along  [1 12]  and  the  c-axis  along  [1  1  1]. 
The  dimensions  of  the  simulation  cell  were  62.0  nm  along 
the  x-axis  and  31.5  nm  along  both  the  y-  and  r-axes,  corre¬ 
sponding  to  a  simulation  cell  of  5,405,160  atoms.  The  par¬ 
tially  cross-slipped  core  structures  obtained  for  the  screw 
(1  /2[T  1  0]  Burgers  vector)-120°  ( 1  /2[T  0 1]  Burgers  vector, 
[Oil]  line  direction)  intersection  in  Rao  et  al.  [6],  as  well 
as  for  Cu,  were  subjected  to  uniform  compressive  or  tensile 
stresses  applied  along  the  [1  1  1]  direction.  Under  such  a 
state  of  stress  there  were  no  resolved  shear  stresses  acting 
along  the  Burgers  vector  direction  for  either  the  screw  or 
the  120°  dislocation.  Also,  there  were  no  applied  Escaig 
stresses  acting  on  the  edge  components  of  the  screw 
character  dislocation  on  the  (1  1  1)  glide  plane.  However, 
there  were  Escaig  stresses  acting  on  the  edge  components 
of  the  screw  dislocation  on  the  cross-slip  (111)  plane,  as 
well  as  the  forest  dislocation  on  the  (111)  plane.  Fixed 
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Fig.  1.  A  schematic  of  the  simulation  cell  used  in  the  atomistic  simulations  described  in  this  manuscript. 
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boundary  conditions  were  applied  along  all  three  directions 
and  energy  minimization  was  performed  using  the  conju¬ 
gate  gradient  technique.  Such  results  were  analyzed  to  ob¬ 
tain  the  energy  differences  between  the  fully  on  glide  plane 
core  structure,  fully  on  cross-slip  plane  core  structure  and 
the  activated  structures  for  this  intersection. 

The  approach  used  to  obtain  activation  energies  for 
transformation  from  one  core  structure  to  another  was  as 
follows.  As  the  stress  was  varied,  the  relative  magnitudes 
of  Escaig  stresses  on  the  glide  plane  and  cross-slip  plane 
changed,  thus  affecting  a  variation  in  the  length  of  the  par¬ 
tially  cross-slipped  region  of  the  screw  character  disloca¬ 
tion.  Note  that  stabilizing  each  configuration  required 
manipulation  of  the  stresses  before,  after  and  at  the  inflec¬ 
tion  point  in  the  energy-distance  coordinate.  By  releasing 
the  stress  from  the  stabilized  state,  the  configuration  was 
allowed  to  “fall”  towards  either  side  of  the  saddle  point 
configuration.  This  process  allowed  the  identification  of 
the  saddle  point  configuration  (within  the  energy  error  of 
the  technique).  Once  identified,  the  absolute  energy  of  this 
configuration  without  the  applied  stress  was  calculated  by 
freezing  a  small  volume  of  atoms  near  the  core  transforma¬ 
tion  region  and  then  minimizing  the  energy  of  all  other 
atoms  in  the  cell.  This  gave  an  upper  bound  for  the  energy 
of  the  saddle  point.  The  activation  energy  was  obtained  by 
subtracting  the  energy  of  the  fully  glide  plane  spread  core 
configuration  from  that  of  the  saddle  point  configuration. 
The  primary  error  in  this  calculation  arose  from  the  fact 
that  an  arbitrary  number  of  atoms  were  kept  fixed  while 
obtaining  the  energy  of  the  saddle  point  configuration 
without  applied  stress  and  was  estimated  to  be  of  the  order 
of  0.2  eV. 

2.1.  Interatomic  potential 

The  embedded  atom  potentials  used  for  the  simulations 
were  those  developed  for  face-centered  cubic  (fee)  Ni  by 
Angelo,  Moody  and  Baskes  [8]  based  on  the  Voter  and 
Chen  format  (Angelo),  as  well  as  the  two  Ni  potentials 
developed  based  on  the  Voter  and  Chen  format  which 
Rao  et  al.  used  in  their  bulk  cross-slip  simulations,  vni 
and  vnih  [9].  The  embedded  atom  potential  used  for  simu¬ 
lations  of  Cu  was  that  developed  by  Mishin  et  al.  [10]. 
Table  1  gives  the  lattice  parameter,  cohesive  energy,  elastic 
constants  and  stacking  fault  energy  for  each  of  the  poten- 

Table  1 


Lattice  parameter  (a0),  cohesive  energy  (Ec),  elastic  constants  (Cu,  C\2 
and  C44)  and  stacking  fault  energy  (7)  given  by  the  three  EAM  potentials 
for  Ni  Angelo,  vni  and  vnih,  as  well  as  the  Mishin  EAM  potential  for  Cu. 


Property 

vni 

vnih 

Angelo 

Mishin 

a0  (nm) 

0.3518 

0.3526 

0.3520 

0.3615 

Cu  (Nm-2) 

2.44  x  1011 

2.44  x  1011 

2.46  x  10" 

1.70  x  10" 

C12  (N  m-2) 

1.49  x  10“ 

1.49  x  1011 

1.47  x  10" 

1.23  x  1011 

C44  (N  rrT2) 

1.25  x  1011 

1.19  x  1011 

1.25  x  10" 

0.76  x  10" 

Ec  (eV) 

-4.43 

-4.43 

-4.45 

-3.54 

y  (J  n-r2) 

00.059 

0.119 

0.089 

0.044 

tials.  The  three  Ni  potentials  used  in  the  simulations  give 
almost  identical  elastic  constants,  cohesive  energies  and  lat¬ 
tice  parameters  (close  to  experimental  values),  whereas  the 
stacking  fault  energy  given  by  the  potentials  varied  from  60 
to  120  mJ  m-2.  The  Shockley  partial  spacing  width  d  for 
the  screw  dislocation  varies  from  d/b  =  5  for  potential  vnih 
to  d/b  =  8  for  potential  vni.  The  Angelo  potential  gives  a  d/ 
b  ratio  of  approximately  six.  This  allows  determination  of 
the  cross-slip  activation  energy  as  a  function  of  d/b  ratio. 
Since  the  experimentally  determined  stacking  fault  energy 
of  Ni  is  close  to  120  mJ  m~2,  the  results  for  cross-slip  acti¬ 
vation  energy  from  potential  vnih  should  be  representative 
of  Ni.  The  EAM  potential  developed  for  Cu  by  Mishin 
et  al.  [10]  gives  lattice  parameter,  cohesive  energy,  elastic 
constants  and  stacking  fault  energy  values  which  are  very 
near  the  experimental  values  and  cross-slip  activation 
energy  results  from  this  potential  should  be  representative 
of  Cu.  The  Mishin  potential  gives  a  d/b  ratio  of  6  for  Cu. 
This  is  identical  to  the  d/b  ratio  given  for  Cu  by  the  poten¬ 
tial  Rasmussen  et  al.  [11]  used  in  their  simulations  of  the 
Escaig  process  for  cross-slip  in  bulk. 

2.2.  Depiction  of  core  structures 

In  order  to  illustrate  the  relaxed  screw  dislocation  geom¬ 
etries  we  took  advantage  of  the  increase  in  atomic  energy 
produced  by  the  strain  field  of  the  partial  dislocations.  By 
plotting  the  atoms  having  assigned  energies  within 
LAAMPs  of  greater  than  —4.42  or  —4.40  eV  (Ni)  or 
-3.52  eV  (Cu)  (the  energy  of  atoms  in  the  stacking  fault 
region)  the  Shockley  partial  dislocation  cores  can  be  imaged 
easily,  even  for  these  large  simulation  cells.  In  order  to  illus¬ 
trate  the  cross-slipped  segment  products  of  the  screw  dislo¬ 
cation  the  positions  are  shown  in  a  [1  1  1]  as  well  as  a  [1 1 2] 
projection  (viewed  directions).  For  the  [1 1 2]  projection  seg¬ 
ments  spread  on  the  initial  (111)  plane  appear  as  a  single 
line  and  cross-slipped  segments  (i.e.  on  a  (111)  plane) 
appear  as  a  pair  of  partials  separated  by  a  stacking  fault. 

3.  Review  of  atomistic  core  structures  for  the  120° 
intersection 

Figs.  2-4  depict  the  [1  1  1]  and  [1 1 2]  projections  (i.e.  the 
x-y  and  x-z  planes)  of  the  various  core  structures  obtained 
for  the  screw  dislocation-1200  dislocation  intersection,  dis¬ 
cussed  in  Rao  et  al.  [6],  obtained  using  the  Angelo  EAM 
potential  for  Ni.  The  starting  configurations  and  boundary 
conditions  were  discussed  in  detail  in  our  previous  paper 
[6].  Fig.  2  illustrates  a  screw  dislocation  fully  on  the  glide 
(111)  plane  and  constricted  (i.e.  Stroh  constrictions  [12]) 
at  both  the  intersections  around  one  of  the  Shockley  par¬ 
tials  of  the  120°  1/2[1 01]  dislocation  residing  on  the 
(111)  plane.  Fig.  3  depicts  two  Lomer-Cottrell  lock  config¬ 
urations  (LC1  and  LC2)  in  the  [1  1  1]  and  [1 1 2]  projections 
obtained  for  the  screw  dislocation- 120°  dislocation  inter¬ 
section.  After  relaxation  the  screw  dislocation  fully  resides 
on  the  cross-slip  (111)  plane  (Fig.  3a  and  b).  In  Fig.  3b,  as 
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GL1 


Fig.  2.  [1  1  1]  (x-y  plane)  and  [1 12]  (x-z  plane)  projections  of  the  core 
structure  (GL1)  for  a  screw-1200  intersection  in  fee  Ni  obtained  from 
atomistic  simulations  using  the  Angelo  EAM  potential.  Atoms  having  an 
energy  greater  than  the  energy  at  a  stacking  fault  are  plotted  and 
correspond  to  the  cores  of  Shockley  partial  dislocations. 

in  Fig.  2,  the  screw  dislocation  is  constricted  at  both  the 
intersections  around  one  of  the  Shockley  partials  of  the 
120°  1/2[1 01]  dislocation  residing  on  the  (111)  plane.  In 
Fig.  3a  the  a/6[l  2 1]  Shockley  partial  of  the  screw  disloca¬ 
tion  combines  with  the  a/6[ll2]  Shockley  partial  of 
the  120°  dislocation  to  form  an  extended  node  resulting 
in  a  short  segment  having  a  Burgers  vector  of  the  type 
l/6[0l  1]  (e.g.  a  Lomer-Cottrell  barrier).  At  the  other  end 
of  the  Lomer-Cottrell  barrier,  the  screw  dislocation  is  con¬ 


stricted  to  form  one  half  of  a  Stroh  constriction.  Fig.  4 
shows  the  [1  1  1]  and  [1 12]  projections  for  two  core  struc¬ 
tures  identified  as  partially  cross-slipped  locks  (PCS1  and 
PCS2)  obtained  for  the  screw  dislocation-1200  dislocation 
intersection.  In  both  cases  the  screw  dislocation  is  partially 
on  the  glide  (111)  plane  and  partially  on  the  cross-slip 
(111)  plane.  In  Fig.  4a,  as  in  Fig.  2,  the  screw  dislocation 
is  constricted  at  both  the  intersections  (negative  constric¬ 
tions  [9,11])  around  one  of  the  Shockley  partials  of  the 
120°  1/2[1 0 1]  dislocation  residing  on  the  (1 1 1)  plane.  Since 
the  core  is  in  the  partially  cross-slipped  state,  a  companion 
positive  constriction  [9,1 1]  forms  between  the  intersections. 
In  Fig.  4b,  as  in  Fig.  3a,  the  a/6[l  2 1]  Shockley  partial  of 
the  screw  dislocation  combines  with  the  a/6[l  1 2]  Shockley 
partial  of  the  120°  dislocation  to  form  an  extended  node 
resulting  in  a  short  segment  having  a  Burgers  vector  of  type 
1/6[01 1]  (e.g.  a  Lomer-Cottrell  barrier).  At  the  other  end 
of  the  Lomer-Cottrell  barrier  the  screw  dislocation  is  con¬ 
stricted  to  form  one  half  of  a  negative  constriction  and  is 
cross-slipped  onto  the  glide  (111)  plane.  Since  the  core  is 
in  the  partially  cross-slipped  state,  a  companion  full  posi¬ 
tive  constriction  forms  between  the  intersections  [9, 1 1  ]. 

Fig.  5  depicts  the  [1  1  1]  and  [1 1 2]  projections  (i.e.  the  x- 
y  and  x-z  planes)  of  two  of  the  core  structures  obtained  for 
the  screw  dislocation-1200  dislocation  intersection, 
obtained  using  the  Mishin  EAM  potential  for  Cu.  The  core 
structure  shown  in  Fig.  5b  for  Cu  is  very  similar  to  the  core 
structure  shown  in  Fig.  2  for  Ni,  with  the  screw  dislocation 
fully  on  the  glide  (111)  plane  GL1  and  the  core  structure 
shown  in  Fig.  5a  for  Cu  is  very  similar  to  the  core  structure 
shown  in  Fig.  4a  for  Ni,  with  the  screw  dislocation  partially 
on  the  glide  (111)  plane  and  partially  on  the  cross-slip 
(111)  plane  PCS1. 


[1-10]  (A) 


[1-10]  (A) 


[1-10]  (A)  [1-10]  (A) 


Fig.  3.  [Ill]  (x-y  plane)  and  [1 1  2]  (x-z  plane)  projections  of  the  core  structures  (LC1  and  LC2)  for  a  screw-1200  intersection  in  fee  Ni  from  atomistic 
simulations  using  the  Angelo  EAM  potential. 
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(a)  pcsi  (b)  PCS2 


Fig.  4.  [1  1  1]  (x-y  plane)  and  [1  1  2]  (x-z  plane)  projections  of  the  core  structures  (PCSI  and  PCS2)  for  a  screw-1200  intersection  in  fee  Ni  from  atomistic 
simulations  using  the  Angelo  EAM  potential. 


x:  [1-10]  (A)  x:  [1-10]  (A) 


x:  [1-10]  (A)  x:  [1-10]  (A) 


Fig.  5.  [1  1  1]  ( x-y  plane)  and  [1  1 2]  (x-z  plane)  projections  of  two  of  the  core  structures  for  a  screw-1200  intersection  in  fee  Cu  from  atomistic  simulations 
using  the  Mishin  EAM  potential. 


3.1.  Schematic  description  of  the  energy  well 

Fig.  6  gives  a  schematic  plot  of  the  energy  of  the  screw 
dislocation  configuration  at  the  intersection  as  a  function 
of  the  length  of  screw  dislocation  cross-slipped  onto  the 
cross-slip  plane  in  between  the  120°  intersecting  disloca¬ 
tions.  Previous  molecular  statics  simulation  results  [6] 
showed  that  there  are  local  minima  in  the  energy  curve 
when:  the  screw  dislocation  length  on  the  cross-slip  plane 
is  zero  (fully  glide  state),  the  entire  screw  dislocation  length 
is  on  the  cross-slip  plane  (fully  cross-slipped  state)  and 


when  approximately  half  the  screw  dislocation  length 
between  the  intersections  is  on  the  cross-slip  plane  (partial 
cross-slip  state).  Two  local  maxima  are  expected  for  the 
energy  versus  state  curve,  in  between  the  fully  glide  and 
partial  cross-slip  states  and  the  fully  cross-slip  and  partial 
cross-slip  states.  Starting  from  the  partial  cross-slip  state, 
application  of  increasing  Escaig  stresses  to  the  configura¬ 
tion  were  used  to  drive  the  screw  dislocation  up  the  energy 
well  until  a  point  of  inflexion  (e.g.  the  first  derivative  of  the 
energy  curve  and  the  stress  was  maximum).  For  stresses 
slightly  above  the  point  of  inflexion  in  the  energy  curve 
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Fig.  6.  (a)  A  schematic  plot  of  the  energy  profile  of  the  screw  dislocation  configuration  as  a  function  of  the  length  of  the  screw  dislocation  on  the  cross-slip 
plane  in  between  the  forest  dislocation  intersections,  (b)  A  schematic  sketch  of  the  [112]  projection  of  the  screw  dislocation  configuration  in  between  forest 
dislocation  intersections  under  applied  Escaig  stresses. 


the  screw  dislocation  will  collapse  to  either  the  fully  glide 
or  fully  cross-slip  state.  The  work  done  by  the  Escaig  stres¬ 
ses  up  to  the  inflexion  point  gives  the  energy  required  for 
the  screw  dislocation  to  transfer  from  the  partially  cross- 
slipped  state  to  the  fully  glide  plane  or  fully  cross-slip  plane 
states.  The  maximum  Escaig  stresses  under  which  the  par¬ 
tial  cross-slipped  state  is  stable  are  given  by  the  maximum 
of  the  first  derivative  of  this  energy  profile. 

The  following  procedure  was  used  to  determine  the  local 
maxima,  i.e.  the  activated  state,  for  transforming  the  screw 
dislocation  from  fully  glide  plane  or  fully  cross-slip  plane 
states  to  the  partially  cross-slipped  state  in  the  energy 
curve.  Partially  cross-slipped  states  under  Escaig  stresses 
slightly  above  the  maximum  of  the  first  derivative  of  the 
energy  curve  were  relaxed  and  stopped  at  different  points 
along  the  energy  curve.  The  Escaig  stresses  were  relaxed 
back  to  zero.  Depending  on  the  initial  value  of  the  Escaig 
stress  and  the  position  on  the  activation  surface,  under  zero 
stress  the  screw  dislocation  configuration  will  collapse  to 
one  of  the  three  possible  minimum  energy  states:  partially 
cross-slipped,  fully  glide  plane  or  fully  cross-slipped.  The 
point  at  which  the  screw  dislocation  transitions  from  relax¬ 
ing  to  the  partially  cross-slipped  state  to  relaxing  to  the 
fully  glide  plane  or  fully  cross-slipped  plane  states  is  iden¬ 
tified  as  the  maxima  in  the  energy  curve  or  the  activated 
state  for  cross-slip.  The  energy  of  this  configuration  relative 
to  the  fully  glide  plane  or  fully  cross-slip  plane  states  was 


evaluated  by  using  two  relaxation  techniques  from  the  acti¬ 
vated  configuration,  which  are  described  in  Section  3.4. 

3.2.  Behavior  of  cores  PCS1  and  PCS2  under  stress 

Fig.  7  gives  the  relaxed  core  structures  in  the  [112]  pro¬ 
jection  for  the  screw-1200  intersection  (core  structure 
PCS1)  for  uniform  compressive  applied  stresses  of  250, 
500,  750.  1000,  1100  and  1150  MPa  along  the  [1  1  1]  direc¬ 
tion,  obtained  using  the  Angelo  EAM  potential.  For  the 
1 1 50  MPa  case  an  additional  increment  of  50  MPa  com¬ 
pressive  stress  was  applied  on  the  relaxed  1 100  MPa  com¬ 
pressive  stress  configuration.  With  increasing  compressive 
applied  stress  the  amount  of  cross-slipped  screw  segment 
in  between  the  two  intersecting  120°  dislocations  monoton- 
ically  increased  until  a  stress  of  1100  MPa  was  reached. 
When  the  compressive  applied  stress  was  slightly  increased 
from  1100  to  1150  MPa  the  screw  dislocation  in  between 
the  intersections  suddenly  collapsed  to  core  structure 
LC2,  fully  residing  on  the  cross-slip  (111)  plane.  For  com¬ 
pressive  applied  stresses  the  effective  stacking  fault  energy 
on  the  cross-slip  (111)  plane  decreased  by  an  amount  equal 
to  t ebe,  where  te  =  Sac,  S  is  the  Schmidt  factor  acting  on 
the  edge  component  of  the  Shockley  partials  on  the 
cross-slip  (111)  plane,  <xc  is  the  applied  compressive  stress 
along  the  [111]  direction  and  be  is  the  magnitude  of  the 
edge  component  =  0.071  nm. 
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(a)  Relaxed 
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(b)  250  MPa  (c)  500  MPa 


x  11-10]  (A)  x  1 1  -1 0]  (A) 


(d)  750  MPa 


x  [1-10]  (A) 


(e)  1000  MPa 
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(f)  1100  MPa 
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Fig.  7.  The  [1  12]  projection  of  the  equilibrium  configurations  of  the  partially  cross-slipped  core  structure,  PCS1,  under  uniform  applied  compressive 
stresses  of  0,  250,  500,  750,  1000,  1 100  and  1150  MPa  applied  along  the  [111]  direction,  obtained  using  the  Angelo  EAM  potential  for  Ni. 


Similar  behavior  was  observed  for  uniform  tensile 
applied  stresses  along  the  [111]  direction.  In  this  case, 
the  amount  of  cross- slipped  segment  in  between  the  two 
intersecting  120°  dislocations  montonically  decreased  until 
a  stress  of  1000  MPa  was  reached.  When  the  tensile  applied 
stress  was  slightly  increased  from  1000  to  1050  MPa  the 
screw  dislocation  in  between  the  intersections  suddenly 
relaxed  to  core  structure  GL1,  fully  residing  on  the  glide 
(1  1  1)  plane.  For  tensile  applied  stresses  the  effective  stack¬ 
ing  fault  energy  on  the  cross-slip  (111)  plane  increased  by 
an  amount  equal  to  t ebe,  where  te  =  Sa,  and  at  is  the 
applied  tensile  stress  along  the  [1  1  1]  direction. 


Similar  behavior  was  observed  for  core  structure  PCS2. 
In  this  case,  when  the  applied  stress  was  compressive  core 
structure  PCS2  relaxed  to  core  structure  LC1  when  the 
compressive  stress  was  increased  from  650  to  700  MPa. 
Similarly,  for  tensile  applied  stresses  core  structure  PCS2 
relaxed  to  core  structure  GL1  when  the  applied  tensile 
stress  was  increased  from  1150  to  1200  MPa. 

3.3.  Analysis  of  cores  PCS1  and  PCS2  under  applied  stresses 

Fig.  8a  and  b  illustrates  the  change  in  length  of  the 
cross-slipped  segment  in  between  the  two  intersecting 
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Fig.  8.  A  plot  of  the  increase  or  decrease  in  length  of  the  cross-slipped 
segment  in  between  the  two  intersecting  120°  dislocations  (x-axis)  versus 
the  decrease  or  increase  in  effective  stacking  fault  energy  (y-axis)  due  to  the 
applied  compressive  or  tensile  stresses  along  the  [111]  direction  for  initial 
core  structure  PCS1  obtained  using  the  Angelo  EAM  potential  for  Ni.  The 
activation  barrier  for  the  transition  of  the  dislocation  from  the  partially- 
cross-slipped  state  (PCS1)  to  the  either  the  fully  cross-slipped  plane  state 
(LC2)  or  the  fully  glide  plane  state  (GL1)  is  also  indicated. 


120°  dislocations  (x-axis)  versus  the  change  in  effective 
stacking  fault  energy  (y-axis)  due  to  the  applied  compres¬ 
sive  or  tensile  stresses  along  the  [111]  direction  for  initial 
core  structure  PCS1  obtained  using  the  Angelo  EAM 
potential.  Fig.  8a  and  b  clearly  shows  that  for  a  certain  crit¬ 
ical  increase  or  decrease  in  stacking  fault  energy,  5y,  the 
core  structure  for  the  intersection  collapsed  to  either 
LC2,  fully  residing  on  the  cross-slip  plane  (compressive 
stress),  or  GL1,  fully  residing  on  the  glide  plane  (tensile 
stress).  The  work  done  by  the  applied  stresses  up  to  the 
critical  point  is  given  by  the  area  under  the  curve  up  to 
the  critical  point  multiplied  by  the  Shockley  partial  spacing 
(estimated  to  be  6b,  where  b  is  the  magnitude  of  the  Bur¬ 
gers  vector  of  the  screw  dislocation  [9])  and  is  equal  to 
the  energy  difference  between  the  activated  states  and  the 
partially  cross-slipped  state  PCS1,  <5 Fact  pcsi-  This  is  deter¬ 
mined  to  be  0.68  eV  for  the  activated  state  on  the  compres¬ 
sive  side  and  1 .44  eV  for  the  activated  state  on  the  tensile 
side.  The  Escaig  stresses  required  to  transform  the  disloca¬ 
tion  from  the  partially  cross-slipped  state  (PCSI)  to  the 
fully  glide  plane  state  (on  the  tensile  side  GL1)  is 
314  MPa  and  to  the  fully  cross-slipped  state  (LC2)  is 
346  MPa. 


Similar  analysis  was  performed  for  core  structure 
PCS2.The  Escaig  stresses  required  to  transform  the  dislo¬ 
cation  from  the  partially  cross-slipped  state  (PCS2)  to  the 
fully  glide  plane  state  (on  the  tensile  side  GL1)  is 
361  MPa  and  to  the  fully  cross-slipped  state  (on  the  com¬ 
pressive  side  LC1)  is  204  MPa. 

3.4.  Determination  of  cross-slip  activation  energy 

To  determine  the  activation  energy  for  cross-slip  from 
the  fully  glide  plane  state  (GL1)  to  the  partially  cross- 
slipped  state  (PCSI)  the  simulation  having  1050  MPa  ten¬ 
sile  stress  applied  to  the  partially  cross-slipped  state  PCSI 
was  stopped  after  varying  amounts  of  relaxation  and  the 
applied  stress  relaxed  back  to  zero.  Conjugate  gradient 
minimization  after  the  stress  was  relaxed  back  to  zero  indi¬ 
cated  whether  the  core  relaxed  to  the  partially  cross-slipped 
state  (PCSI)  or  the  fully  glide  plane  state  (GL1).  The  con¬ 
figuration  at  which  the  relaxation  transitions  from  relaxing 
to  PCSI  to  relaxing  to  GL1,  was  taken  to  be  the  activated 
state  for  cross-slip  from  GL1  to  PCSI.  Fig.  9  gives  the 
structure  of  the  activated  state  in  the  [1 1 2]  projection  for 
the  screw-1200  intersection  (core  structure  PCSI)  obtained 


(a)  Tension 


(b)  Compression 


Fig.  9.  The  [1  1 2]  projection  of  the  structure  of  the  activated  configuration 
for  transition  of  the  screw  dislocation  from  (a)  state  GL1  to  PCSI  (tensile 
side)  and  (b)  state  LC2  to  state  PCSI  (compressive  side)  obtained  using  the 
Angelo  EAM  potential  for  Ni.  The  x-axis  is  along  the  [1  1 0]  direction  and 
the  j>-axis  is  along  the  [111]  direction.  The  square  boxes  indicate  regions 
shown  at  higher  magnification  in  the  insets. 
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using  the  Angelo  potential,  which  gives  a  stacking  fault 
energy  of  90  mJ  m“  2.  Fig.  9a  shows  that  in  the  activated 
configuration  a  small  portion  of  the  screw  dislocation  near 
the  intersection,  of  length  ~1.5<i,  resides  on  the  cross-slip 
plane.  The  attractive  interaction  between  the  negative 
and  the  positive  constrictions  becomes  zero  at  a  separation 
distance  of  1.5 d,  beyond  which  it  becomes  repulsive,  in  the 
presence  of  an  intersecting  dislocation.  This  behavior  is  sig¬ 
nificantly  different  from  the  typical  Escaig  process  for 
cross-slip,  where  the  attractive  interaction  between  the  con¬ 
strictions  only  becomes  zero  at  an  infinite  separation  dis¬ 
tance  [2-4]. 

A  similar  analysis  was  performed  to  determine  the  acti¬ 
vation  energy  for  cross-slip  from  the  fully  cross-slip  plane 
state  (LC2)  to  the  partially  cross-slipped  state  (PCS1).  In 
this  case  the  simulation  having  a  1150  MPa  compressive 
stress  applied  to  the  partially  cross-slipped  state  PCS1 
was  stopped  after  varying  amounts  of  relaxation  and  the 
applied  stress  relaxed  back  to  zero.  Conjugate  gradient 
minimization  after  the  stress  was  relaxed  back  to  zero  indi¬ 
cated  whether  the  core  relaxed  to  the  partially  cross-slipped 
state  (PCS1)  or  the  fully  cross-slip  plane  state  (LC2).  The 
configuration  at  which  the  relaxation  transitioned  from 
relaxing  to  PCS1  to  collapsing  to  LC2  was  taken  to  be 
the  activated  state  for  cross-slip  from  LC2  to  PCS1. 
Fig.  9b  shows  that  in  the  activated  configuration  a  small 
portion  of  the  screw  dislocation  near  the  intersection,  of 
length  resided  on  the  glide  plane.  The  attractive  inter¬ 
action  between  the  negative  and  the  positive  constrictions 
became  zero  at  a  separation  distance  of  d,  beyond  which 
it  became  repulsive,  in  the  presence  of  an  intersecting  dislo¬ 
cation.  Again,  such  behavior  is  significantly  different  from 
the  standard  Escaig  process. 

A  similar  analysis  was  performed  on  the  partially  cross- 
slipped  state  (PCSlj  obtained  using  the  Mishin  potential 
for  Cu,  and  the  activated  states  obtained  using  this  poten¬ 
tial  on  both  the  tensile  and  compressive  sides  were  very 
similar  to  the  activated  configurations  obtained  using  the 
Angelo  potential  for  Ni. 

The  activation  energy  for  cross-slip  from  state  GL1  to 
state  PCS1,  b/f(,M  pc  s i  was  determined  using  two  different 
relaxation  techniques  from  the  identified  activated  configu¬ 
ration.  In  the  first  technique  the  energy  was  determined  for 
relaxation  from  the  activated  state  to  state  GL1  under  an 
applied  tensile  stress  of  1050  MPa  (Angelo  potential  for 
Ni)  along  the  [1  1  1]  direction.  This  energy  was  corrected 
for  two  factors  to  obtain  the  activation  energy  for  cross¬ 
slip:  (a)  work  done  by  the  applied  Escaig  stresses  in  collaps¬ 
ing  the  activated  configuration  (of  the  order  of  0.2  eV);  (b) 
the  increase  in  energy  of  the  activated  configuration  due  the 
applied  Escaig  stresses  (of  the  order  of  0.4  eV).  For  the  sec¬ 
ond  technique  the  activated  state  was  relaxed  to  state  GL1 
having  zero  applied  stress,  with  and  without  atoms  within  a 
cube  of  size  1.2  nm  centered  at  the  positive  constriction  of 
the  activated  state  fixed.  The  difference  in  relaxation  ener¬ 
gies  between  the  two  types  of  simulations  was  taken  to 
be  the  activation  energy  for  cross-slip  from  GL1  to 


PCS1.  The  results  for  the  cross-slip  activation  energy 
obtained  using  these  two  techniques  of  relaxation  were 
within  0.2  eV  of  each  other  and  were  taken  to  be  represen¬ 
tative  of  the  error  in  the  determination  of  the  activation 
energy  for  cross-slip  using  this  procedure.  The  activation 
energy  for  cross-slip  from  state  GL1  to  state  PCS1, 
5EGU  pcs i ,  was  also  determined  when  the  separation  dis¬ 
tance  between  the  two  intersecting  forest  dislocations  was 
90  nm,  as  opposed  to  30  nm,  and  the  cross-slip  activation 
energy  was  found  to  be  independent  of  the  separation  dis¬ 
tance  used  in  the  simulations. 

Similar  simulations  and  analyses  were  performed  to 
determine  the  activation  energy  for  cross-slip  from  GL1 
to  PCS1  using  two  other  Ni  EAM  potentials,  vni  and  vnih. 
Also,  similar  simulations  and  analyses  were  performed  to 
determine  the  activation  energy  for  cross-slip  from  the  fully 
cross-slip  plane  state  (LC2)  to  the  partially  cross-slipped 
state  (PCS  1 )  using  all  three  EAM  potentials,  dElC2  pc  si- 
For  these  cases  the  two  techniques  of  relaxation  described 
before  for  determining  the  cross-slip  activation  energy  gave 
values  that  differed  by  less  than  0.1  eV. 

Fig.  10  is  a  plot  of  the  cross-slip  activation  energies  for 
cross-slip  from  state  GL1  to  state  PCSI  as  well  as  from 
state  LC2  to  state  PCSI,  obtained  using  the  three  Ni 
EAM  potentials,  versus  the  parameter  (r//b)[ln(v/3r//b)]1/2 
given  by  the  potentials.  Also  plotted  in  Fig.  10  are  the  cor¬ 
responding  values  for  the  cross-slip  activation  energy  deter¬ 
mined  from  simulations  of  the  Escaig  process  as  well  as  the 
cross-slip  activation  energies  obtained  for  the  case  of  a  sin¬ 
gle  intersecting  dislocation,  as  opposed  to  two,  with  the 
Angelo  Ni  potential  [9].  Fig.  10  clearly  shows  that  the 
cross-slip  activation  energy  is  linearly  proportional  to  ( d / 
b)[ln(  v/3r//b)]1/2  both  at  the  intersection  and  for  the  Escaig 
type  process.  The  cross-slip  activation  energy  (from  GL1  to 
PCSI)  values  near  the  intersection  were  a  factor  of  three 
lower  than  that  for  the  Escaig  process.  Also,  the  cross-slip 


(d/b)[ln(31,2d/b)]1/2 

Fig.  10.  A  plot  of  the  activation  energy  for  cross-slip  from  state  GL1  to 
state  PCSI  (^-Egli-pcsi)  as  well  as  from  state  LC2  to  state  PCSI  (SELci- 
pcsi)  in  Ni  as  a  function  of  the  parameter  J/b[ln(\/3r//b)]l/2  given  by  the 
EAM  potentials  Angelo,  vni  and  vnih.  Also  shown  are  the  corresponding 
cross-slip  activation  energy  values  for  the  Escaig  process  operating  in  a 
bulk  crystal  [9]  as  well  as  cross-slip  activation  energy  values  obtained  for 
the  case  of  a  single  intersecting  dislocation,  as  opposed  to  two,  with  the 
Angelo  potential. 
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activation  energy  from  the  fully  cross-slip  plane  state  (LC2) 
to  the  partially  cross-slipped  state  (PCS1)  was  approxi¬ 
mately  a  factor  of  five  lower  than  that  for  the  Escaig  pro¬ 
cess.  The  <5ELC2_pcsi  value  can  be  identified,  by 
symmetry,  as  the  cross-slip  activation  energy  from  the  fully 
glide  plane  state  to  the  partially  cross-slipped  state  at  a 
120°  screw  dislocation  intersection,  where  the  intersecting 
forest  dislocation  has  a  1/2(0  1  1)  Burgers  vector  ( 1 0 1), 
line  direction  and  resides  on  the  (ill)  plane.  This  suggests 
that  the  cross-slip  activation  energy  near  an  intersection  is 
strongly  dependent  on  the  Burgers  vector  of  the  intersect¬ 
ing  dislocation  and  could  be  a  factor  of  five  lower  than  that 
usually  considered.  Fig.  10  also  shows  that  the  presence  of 
a  second  intersecting  dislocation  has  very  little  effect  on  the 
calculated  cross-slip  activation  energies  at  the  intersections. 

The  activation  volume  for  the  cross-slip  process  with 
respect  to  applied  Escaig  stresses  is  given  by  the  area  of 
the  cross-slipped  portion  of  the  screw  character  dislocation 
in  the  activated  state  (Fig.  9)  multiplied  by  the  magnitude 
of  the  edge  component  of  the  Burgers  vector  of  the  Shock- 
ley  partials.  From  the  simulations  values  for  this  parameter 
ranged  from  10  to  20b3  for  cross-slip  at  screw  dislocation 
intersections. 

Similarly,  the  cross-slip  activation  energies  at  the  120° 
intersection,  dFom-pcsi  as  well  as  dE\<  2  pc  si-  were  deter¬ 
mined  using  the  Mishin  potential  for  Cu.  These  values  were 
determined  to  be  1.13  and  0.65  eV,  respectively.  These  val¬ 
ues  are  factors  of  2-3  lower  than  that  for  the  Escaig  process 
(1.7  eV)  (unpublished  calculations  of  cross-slip  activation 
energy  via  the  Escaig  process  in  bulk  for  Cu  using  the  Mis¬ 
hin  potential  by  SI  Rao  (positive  constriction  3.0  eV,  neg¬ 
ative  constriction  -1.3  eV)).  The  values  of  1.13  and 
0.65  eV  determined  for  the  cross-slip  activation  energies 
at  the  120°  intersection  are  in  reasonable  accord  with  the 
experimentally  determined  value  for  cross-slip  activation 
energy  in  Cu,  1.15  eV  [13].  This  suggests  that  cross-slip 
preferentially  occurs  at  selected  screw  dislocation  intersec¬ 
tions  in  fee  materials. 

4.  Discussion 

This  work  has  shown  that  the  activation  energy  for 
nucleation  of  cross-slip  in  the  bulk  of  an  fee  metal  can  be 
sufficiently  small  to  enable  profuse  cross-slip  through  ther¬ 
mal  activation  under  multi-slip  conditions.  Further,  this 
work  provides  a  methodology  to  obtain  the  variation  in 
nucleation  barrier  with  type  of  slip  system  interaction 
within  a  given  system.  These  calculations  can  be  directly 
used  in  higher  level  mesoscale  simulations  where  cross-slip 
statistics  can  be  accounted  for.  The  results  are  likely  to  per¬ 
mit  improved  simulations  of  metallic  behavior,  such  as 
work  hardening  and  fatigue,  that  have  not  been  possible 
without  ad  hoc  assumptions  about  cross-slip.  We  include 
below  a  few  remarks  on  our  results,  including  limitations. 

The  results  clearly  show  that  cross-slip  nucleation  at 
screw  dislocations  intersecting  forest  dislocations  exhibits 
an  activation  energy  that  is  a  factor  of  2-5  lower  than  that 


for  the  Escaig  mechanism.  It  is  worth  noting  that  the  prob¬ 
ability  of  dislocation  intersections  is  relatively  large  as  com¬ 
pared  with  other  mechanisms  of  obstacle-induced  cross¬ 
slip,  like  cross-slip  at  screw  dipoles,  jogs  and  surfaces 
[14-16].  Further,  cross-slip  at  screw  dipoles,  jogs  and  sur¬ 
faces  is  athermal  only  under  specific  conditions  [14—16]. 
The  current  results  rationalize  the  experimentally  observed 
profuse  nature  of  cross-slip  in  fee  crystals  in  a  more  satis¬ 
factory  manner  than  the  current  models  of  thermally  acti¬ 
vated  cross-slip  that  require  high  stress  to  provide  a  self- 
consistent  explanation. 

Acknowledging  a  limitation  of  this  work,  we  note  that 
since  a  molecular  statics  technique  has  been  used,  the 
cross-slip  activation  energy  values  reported  in  this  manu¬ 
script  are  expected  to  be  upper  bounds.  A  more  sophisti¬ 
cated  molecular  dynamics  reaction  pathway  technique 
[17]  might  provide  more  accurate  values  for  the  activation 
energies  of  cross-slip  near  forest  dislocation  intersections. 

The  new  cross-slip  nucleation  model  has  a  variety  of 
implications  for  crystal  plasticity  in  fee  materials.  For  exam¬ 
ple,  within  the  present  mechanism  the  frequency  of  cross-slip 
should  scale  with  the  forest  dislocation  density,  i.e.  multipli¬ 
cation  rate  proportional  to  p{.  The  growth  of  such  nuclei 
should  depend  upon  the  relative  magnitude  of  local  stresses 
on  the  glide  plane  and  the  cross-slip  plane  at  the  partially 
cross-slipped  screw-dislocation  intersection  region.  Such 
behavior  of  the  partially  cross-slipped  core  under  different 
modes  of  applied  stress  could  be  studied  using  atomistic  sim¬ 
ulations  as  well  as  dislocation  dynamics  simulations.  Similar 
activation  analyses  using  atomistic  simulations  must  be  per¬ 
formed  for  other  intersections  (varying  the  line  direction  of 
the  intersecting  dislocation)  on  the  (111)  plane  that  form 
glide,  Fomer-Cottrell  or  Hirth  locks  to  determine  the  depen¬ 
dence  of  the  activation  energy  for  cross-slip  on  the  line  direc¬ 
tion  as  well  as  the  Burgers  vector  of  the  intersecting  forest 
dislocation.  The  intersection  mechanism  of  cross-slip  nucle¬ 
ation  should  also  be  implemented  in  3D  dislocation  dynam¬ 
ics  simulations  as  an  alternative  to  Escaig’s  model  for  fee 
materials.  Finally,  it  would  be  instructive  to  revisit  the  Bon¬ 
neville  and  Escaig  experimental  results  on  cross-slip  in  fee  Cu 
in  the  light  of  the  new  intersection  mechanism  for  cross-slip 
nucleation  [1,13]. 

According  to  Washburn  [18],  double  intersection  cross¬ 
slip,  where  the  segment  that  has  been  pulled  into  the  cross¬ 
slip  plane  soon  encounters  another  attractive  intersection 
that  brings  it  back  onto  another  primary  glide  plane,  pro¬ 
vides  a  reasonable  mechanism  for  dislocation  multiplica¬ 
tion  and  the  growth  of  slip  bands  at  low  temperatures. 
Naturally,  such  a  mechanistic  process  must  be  demon¬ 
strated  via  modern  simulation  methods.  Also,  classical  the¬ 
ories  of  strain  hardening  assume  that  dislocation  storage  in 
Stage  II  of  single-crystal  fee  materials  is  a  result  of  junction 
formation  [19]  or  two-dimensional  (2D)  concave  loop  for¬ 
mation  [20]  as  the  gliding  dislocation  traverses  through  an 
array  of  forest  dislocation  obstacles  on  its  glide  plane. 
However,  one  of  the  major  problems  in  classical  strain 
hardening  models  is  to  explain  how  the  generation  of  a 
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3D  network  of  stored  dislocations  occurs  as  a  consequence 
of  2D  glide  [20],  We  note  that  the  intersection  cross-slip 
nucleation  mechanism  for  dislocation  storage  may  provide 
a  convenient  mechanism  of  generating  a  3D  network  of 
stored  dislocations  from  2D  glide. 

Finally,  the  anomalous  flow  behavior  of  intermetallics 
(i.e.  Ni3Al)  has  thus  far  been  modeled  using  core  transfor¬ 
mations  that  require  thermal  activation.  Here  again,  the 
key  issue  has  been  that  constriction  energies  are  prohibi¬ 
tively  high  [2,3,21].  The  possibility  of  core  transformation 
in  the  presence  of  interaction  with  other  dislocations  pro¬ 
vides  a  novel  mechanism  to  consider  and  explore  in  ratio¬ 
nalizing  the  flow  anomaly  of  these  types  of  materials. 

5.  Summary 

Atomistic  simulations  were  conducted  to  identify  the 
saddle  point  configuration  and  thus  the  activation  energy 
for  cross-slip  of  a  screw  dislocation  in  the  presence  of  a  pair 
of  forest  dislocations.  The  approach  involved  the  use  of 
external  stress  along  [1  1  1]  to  vary  the  segment  length  for 
a  partially  cross-slipped  1/2[1 1 0]  screw  dislocation  inter¬ 
secting  a  pair  of  120°  1/2[1 01]  Burgers  vector,  (Oil)  line 
direction  dislocations  residing  on  the  (111)  plane.  The  sim¬ 
ulations  in  this  work  show  that  the  nucleation  of  cross-slip 
at  dislocation-dislocation  interactions  is  significantly  more 
probable  than  that  for  the  Escaig  process,  having  an  activa¬ 
tion  barrier  that  is  a  factor  of  2-5  lower.  Further,  this  work 
shows  that  atomistic  simulations  can  be  used  to  identify 
types  of  interactions  that  result  in  cross-slip  and  these 
results  in  turn  may  be  used  in  higher  level  mesoscale 
simulations. 

Summarizing  the  results,  equilibrium  partially  cross- 
slipped  core  structures  under  applied  compressive  or  tensile 
stresses  along  the  [111]  direction  were  determined  using 
atomistic  simulations.  Beyond  a  critical  stress  the  equilib¬ 
rium  configurations  relax  to  a  core  structure  that  is  either 
fully  on  the  glide  plane  for  tensile  stresses  or  fully  on  the 
cross-slip  plane  for  compressive  stresses. 

The  activation  barrier  for  cross-slip  from  the  fully  glide 
plane  state  (GL1)  to  the  partially  cross-slipped  state 
(PCS1),  (5£gL1-pcsi,  at  these  intersections  was  determined 
to  be  a  factor  of  three  lower  than  that  for  the  conventional 
Escaig-based  mechanism.  While  the  new  cross-slip  mecha¬ 
nism  has  a  lower  activation  energy  it  has  the  same  linear 
relationship  to  (r//b)[ln(v/3r//b)]I/2  as  the  Escaig  process. 

The  activation  barrier  for  cross-slip  from  the  fully  cross¬ 
slip  plane  state  (LC2)  to  the  partially  cross-slipped  state 
(PCS1),  hfsLC2,  at  these  intersections  was  determined  to 
be  approximately  a  factor  of  five  lower  than  for  the  Escaig 
process.  These  b£,LC2-pcsi  values  were  identified,  by  sym¬ 
metry,  as  the  activation  barrier  for  cross-slip  from  the  fully 
glide  plane  state  to  the  partially  cross-slipped  state  at  a 


120°  screw  dislocation  intersection,  where  the  intersecting 
dislocation  has  a  1/2(0  1  1)  Burgers  vector  (10 1),  line 
direction  and  resides  on  the  ( 1 1 1 )  plane. 

Cross-slip  activation  energy  values  near  forest  disloca¬ 
tion  intersections  in  fee  Cu  were  determined  to  be  1.13 
and  0.65  eV.  Though  anticipated  to  be  upper  bound  values, 
these  values  are  in  reasonable  accord  with  the  experimen¬ 
tally  determined  value  for  cross-slip  activation  energy  in 
Cu,  1.15  eV  [13].  This  suggests  that  cross-slip  preferentially 
occurs  at  selected  screw  dislocation  intersections  in  fee 
materials. 
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